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Abstract: We examine the minimal supergravity (mSUGRA) model under the assump- 
tion that the strong CP problem is solved by the Peccei-Quinn mechanism. In this case, 
the relic dark matter (DM) abundance consists of three components: i). cold axions, ii). 
warm axinos from neutralino decay, and in), cold or warm thermally produced axinos. To 
sustain a high enough re-heat temperature (Tr ~ 10^ GeV) for many baryogenesis mech- 
anisms to function, we find that the bulk of DM should consist of cold axions, while the 
admixture of cold and warm axinos should be rather slight, with a very light axino of mass 
~ 100 keV. For mSUGRA with mainly axion cold DM (CDM), the most DM-preferred 
parameter space regions are precisely those which are least preferred in the case of neu- 
tralino DM. Thus, rather different SUSY signatures are expected at the LHC in the case 
of mSUGRA with mainly axion CDM, as compared to mSUGRA with neutralino CDM. 
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1. Introduction 



The cosmic abundance of cold dark matter (CDM) has been recently measured to high 
precision by the WMAP collaboration|jl|, which lately finds 

ncDMh"^ = 0.110 ±0.006, (1.1) 

where = p/pcis the dark matter density relative to the closure density, and h is the scaled 
Hubble constant. No particle present in the Standard Model (SM) of particle physics has 
exactly the right properties to constitute CDM. However, CDM does emerge naturally from 
two compelling solutions to longstanding problems in particle physics. 

The first problem is the strong CP problem [^, for which an elegant solution was 
proposed by Peccei and Quinn many years agoQ, and which naturally predicts the existence 
of a new particle[Q]: the axion o. The axion turns out to be an excellent candidate particle 
for CDM in the universe Q. 

The second problem- the gauge hierarchy problem- arises due to quadratic divergences 
in the scalar sector of the SM. The quadratic divergences lead to scalar masses blowing 
up to the highest scale in the theory {e.g. in grand unified theories (GUTS), the GUT 
scale Mqut — 2 x 10^^ GeV), unless exquisite fine-tuning of parameters is invoked. The 
gauge hierarchy problem is naturally solved by introducing supersymmetry (SUSY) into 
the theory. By including softly broken SUSY, quadratic divergences cancel between fermion 
and boson loops, and only log divergences remain. The log divergence is soft enough that 
vastly different scales remain stable within a single effective theory. In SUSY theories, the 
lightest neutralino emerges as an excellent WIMP CDM candidate. The gravitino of SUSY 
theories is also a good super- WIMP CDM candidate[^. Gravity-mediated SUSY breaking 
models include gravitinos with weak-scale masses. These models experience tension due 
to possible overproduction of gravitinos in the early universe. In addition, late decaying 
gravitinos may disrupt calculations of light element abundances produced by Big Bang 
nucleosynthesis (BBN). This tension is known as the gravitino problem. 

Of course, it is highly desirable to simultaneously account for both the strong CP 
problem and the gauge hierarchy problem. In this case, it is useful to invoke supersymmetric 
models which include the PQ solution to the strong CP problemQ. In a SUSY context, 
the axion field is just one element of an axion supermultiplet. The axion supermultiplet 
contains a complex scalar field, whose real part is the -R-parity even saxion field s(x), and 
whose imaginary part is the axion field a(x). The supermultiplet also contains an i?-parity 
odd spin-^ Majorana field, the axino aQ. The saxion, while being an i?-parity even field, 
nonethless receives a SUSY breaking mass likely of order the weak scale. The axion mass is 
constrained by cosmology and astrophysics to lie in a favored range 10-2 eV> ma > 10"^ 
eV. The axino mass is very model dependent!^, and is expected to lie in the general range 
of keV to GeV. An axino in this mass range would likely serve as the lightest SUSY particle 
(LSP), and is also a good candidate particle for cold dark matter []lC|]. 

In this paper, we investigate supersymmetric models wherein the PQ solution to the 
strong CP problem is also invoked. For definiteness, we will restrict ourselves to examining 
the paradigm minimal supergravity (mSUGRA or CMSSM) model[ll|. We will restrict our 
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work to cases where the hghtest neutrahno x? is the next-to-hghtest SUSY particle (NLSP); 
the case with a stau NLSP has recently been examined in Ref. [^]. Related previous work 
on axino DM in mSUGRA can be found in Ref. 

We will be guided in our analysis also by considering the possibility of including a 
viable mechanism for baryogenesis in the early universe. In order to do so, we will need to 
allow for re-heat temperatures after the inflationary epoch to reach values Tr ^ 10^ GeV. 
We will find that in order to sustain such high re-heat temperatures, as well as generating 
predominantly cold dark matter, we will be pushed into mSUGRA parameter space regions 
that are very different from those allowed by the case of thermally produced neutralino 
dark matter. In addition, we find that very high values of the PQ breaking scale fa/N of 
order 10^^ - 10^^ Q^y are needed, leading to the mSUGRA model with mainly axiom cold 
dark matter, but also with a small admixture of thermally produced axinos, and an even 
smaller component of warm axino dark matter arising from neutralino decays. The favored 
axino mass value is of order 100 keV. We note here recent work on models with dominant 
axion CDM explore the possibility that axions form a cosmic Bose-Einstein condensate, 
which can allow for the solution of several problems associated with large scale structure 
and the cosmic background radiation 1 15]. 

The remainder of this paper is organized as follows. In Sec. we first discuss the 
gravitino problem, and then examine several possible baryogenesis mechanisms: thermal 
and non-thermal leptogenesis and Affleck-Dine leptogenesis. We then examine production 
mechanisms for axion and thermally and non-thermally produced axino dark matter. In 
Sec. ^, we confront the mSUGRA model with the possibility of mixed axion and axino 
cold and warm dark matter. We plot out contours of re-heat temperature Tr, and find 
that parameter space regions with large enough Tr to sustain at least non-thermal lep- 
togenesis favor a sparticle mass spectrum which is actually most disfavored by mSUGRA 
with neutralino cold dark matter. Likewise, the regions of mSUGRA space most favored 
by neutralino CDM are least favored by mixed axion/axino dark matter. This has a large 
impact on the sort of SUSY signatures to be expected at LHC. The requirement of mainly 
axion CDM with Tr ~ 10^ GeV favors rather heavy squarks and sleptons. Thus, we expect 
in this case that LHC signatures will be dominated by gluino pair production followed by 
3-body gluino decays to charginos and neutralinos. In Sec. we present a summary and 
conclusions. 



2. The gravitino problem, leptogenesis and mixed axion/axino dark mat- 
ter 



We adopt the mSUGRA model[ll| as a template model for examining the role of mixed 
axion/axino dark matter in gravity-mediated SUSY breaking models. The mSUGRA pa- 
rameter space is given by 

mo, mi/2, ^0, tan/3, sign{fi), (2.1) 

where mo is the unified soft SUSY breaking (SSB) scalar mass at the GUT scale, mi/2 
is the unified gaugino mass at Mqut, is the unified trilinear SSB term at Mqut and 
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tan/3 = Vu/vd is the ratio of Higgs field vevs at tlie weak scale. The GUT scale gauge and 
Yukawa couplings, and the SSB terms are evolved using renormalization group equations 
(RGEs) from Mqut to 

iTT'weakj &t which point electroweak symmetry is broken radiatively, 
owing to the large top quark Yukawa coupling. At m^gafc) the various sparticle and Higgs 
boson mass matrices are diagonalized to find the physical sparticle and Higgs boson masses. 
The magnitude, but not the sign, of the superpotential fi parameter is determined by the 
EWSB minimization conditions. We adopt the Isasugra subprogram of Isajet for spectra 



generation [16 1 



2.1 Gravitino problem 

In supergravity models, supersymmetry is broken via the superHiggs mechanism. The 
common scenario is to postulate the existence of a hidden sector which is uncoupled to the 
MSSM sector except via gravity. The superpotential of the hidden sector is chosen such 
that supergravity is broken, which causes the gravitino (which serves as the gauge particle 
for the superHiggs mechanism) to develop a mass in^/2 ~ n^-'^/Mpi ~ myj^ak- Here, m is a 
hidden sector parameter assumed to be of order 10^^ GeV. ^ In addition to a mass for the 
gravitino, SSB masses of order m^^ak are generated for all scalar, gaugino, trilinear and 
bilinear SSB terms. Here, we will assume that 7713/2 is larger than the lightest MSSM mass 
eigenstate, so that the gravitino essentially decouples from all collider phenomenology. 

In all SUGRA scenarios, a potential problem arises for weak-scale gravitinos: the 
gravitino problem. In this case, gravitinos G can be produced thermally in the early 
universe (even though the gravitinos are too weakly coupled to be in thermal equilibrium) 
at a rate which depends on the re-heat temperature Tr of the universe. The produced 
G can then decay to various sparticle-particle combinations, with a long lifetime of order 
1 — 10^ sec (due to the Planck suppressed gravitino coupling constant). The late gravitino 
decays occur during or after BBN, and their energy injection into the cosmic soup threatens 
to destroy the successful BBN predictions of the light element abundances. The precise 



constraints of BBN on the gravitino mass and Tr are presented recently in Ref. [18|. One 
way to avoid the gravitino problem in the case where ~ ^ TeV is to maintain a value of 
Tr ^ 10^ GeV. Such a low value of Tr rules out many attractive baryogenesis mechanisms, 
and so here instead we assume that 171^/2 ^ 5 TeV. In this case, the G is so heavy that 
its lifetime is of order 1 sec or less, and the G decays near the onset of BBN. In this case, 
values of Tr as large as 10^ GeV are allowed. 

In the simplest SUGRA models, one typically finds ttiq = m^/2- For more general 



SUGRA models, the scalar masses are in general non-degenerate and only of order 771,3/2 [19 1- 
Here for simplicity, we will assume degeneracy of scalar masses, but with niQ <^ rn^/2- 

2.2 Leptogenesis 

One possible baryogenesis mechanism that elec- 
troweak baryogenesis. However, calculations of successful electroweak baryogenesis within 
the MSSM context seem to require sparticle mass spectra with rrih ~ 120 GeV, and 



^In Ref. Irrfl, a link is suggested between hidden sector parameters and the PQ breaking scale fa- 
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m^^ ~ 125 GeV [pO| . The latter requirement is difficult (though not impossible) to achieve 



in the MSSM, and is also partially excluded by collider searches for light top squarks|21|. 
We will not consider this possibility further. 

An alternative attractive mechanism- especially in light of recent evidence for neutrino 
mass- is thermal leptogenesis [^] . In this scenario, heavy right-handed neutrino states Ni 
(i = 1 — 3) decay asymmetrically to leptons versus anti-leptons in the early universe. The 
lepton-antilepton asymmetry is converted to a baryon-antibaryon asymmetry via sphaleron 
effects. The measured baryon abundance can be achieved provided the re- heat temperature 



Tji exceeds ~ 10^ GeV|24]. The high Tr value needed here apparently puts this mechanism 
into conffict with the gravitino problem in SUGRA theories. 

A related leptogenesis mechanism called non-thermal leptogenesis invokes an alter- 
native to thermal production of heavy neutrinos in the early universe. In non-thermal 
leptogenesis [p^] , it is possible to have lower reheat temperatures, since the Ni may be 
generated via inflaton decay. The Boltzmann equations for the B — L asymmetry have 



been solved numerically in Ref . [ 26 1 . The B — L asymmetry is then converted to a baryon 
asymmetry via sphaleron effects as usual. The baryon-to-entropy ratio is calculated in | p6| , 
where it is found 

»-"'<(lo^)(^)(a^)*- <-) 

where is the inflaton mass and Jg// is an effective CP violating phase which may be of 
order 1. Comparing calculation with data (the measured value of hb/s ~ 0.9 x lO"^'^), a 
lower bound Tr ~ 10'= GeV may be inferred for viable non-thermal leptogenesis via inflaton 
decay. 



A fourth mechanism for baryogenesis is Affleck-Dine [^] leptogenesis [28]. In this ap- 
proach, a flat direction cj)i = (2ff£j)^/^ is identified in the scalar potential, which may have 
a large field value in the early universe. When the expansion rate becomes comparable 
to the SSB terms, the field oscillates, and since the field carries lepton number, coherent 
oscillations about the potential minimum will develop a lepton number asymmetry. The 
lepton number asymmetry is then converted to a baryon number asymmetry by sphalerons 
as usual. Detailed calculations [p^] find that the baryon-to-entropy ratio is given by 

UB 1 \{H)\''Tr 



s 23 rUyMli 



(2.3) 



where {H) is the Higgs field vev, rriy is the mass of the lightest neutrino and Mpi is the 
Planck scale. To obtain the observed value of n^/s, values of Tr ~ 10^ — 10^ are allowed 
for ~ 10-9 - lO"'^ eV. 

Thus, to maintain accord with either non-thermal or Affieck-Dine leptogenesis, along 
with constraints from the gravitino problem, we will aim for axion/axino DM scenarios 
with Tr ~ 10^ - 10^ GeV. 

2.3 Mixed axion/axino dark matter 
2.3.1 Relic axions 

Axions can be produced via various mechanisms in the early universe. Since their life- 
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time (they decay via a — > 77) turns out to be longer than the age of the universe, they 
can be a good candidate for dark matter. Since we will be concerned here with re-heat 
temperatures Tr ~ 10^ GeV < fa/N (to avoid overproducing gravitinos in the early uni- 
verse), the axion production mechanism relevant for us here is just one: production via 
vacuum mis-alignment|^]. In this mechanism, the axion field 0(2;) can have any value 
~ fa at temperatures T ^ Aqcd- As the temperature of the universe drops, the poten- 
tial turns on, and the axion field oscillates and settles to its minimum at —6fa/N (where 
6 = 9 + arg(det mq), 6 is the fundamental strong CP violating Lagrangian parameter 
and rriq is the quark mass matrix). The difference in axion field before and after potential 
turn-on corresponds to the vacuum mis-alignment: it produces an axion number density 

na{t) ~ ^ma{t){aHt)), (2.4) 

where t is the time near the QCD phase transition. Relating the number density to the 
entropy density allows one to determine the axion relic density today |^]: 

^ ,2 1 /6 X 10-6 eV\^/^ , , 



An error estimate of the axion relic density from vacuum mis-alignment is plus-or-minus a 
factor of three. Axions produced via vacuum mis-alignment would constititute cold dark 
matter. However, in the event that (a-^(t)) is inadvertently small, then much lower values 
of axion relic density could be allowed. Additional entropy production at t > tgcD can 
also lower the axion relic abundance. Taking the value of Eq. ( |2.5| ) literally, and comparing 
to the WMAP5 measured abundance of CDM in the universe, one gets an upper bound 
fa/N ~ 5 X 10^^ GeV, or a lower bound ~ 10~^ eV. If we take the axion relic density a 
factor of three lower, then the bounds change to fa/N ^ 1.2 x 10^^ GeV, and rua ~ 4 x 10"^ 
eV. 

2.3.2 Axinos from neutralino decay 

If the d is the lightest SUSY particle, then the x? will no longer be stable, and can decay via 
Xi dj. The relic abundance of axinos from neutralino decay (non-thermal production, 
or NTP) is given simply by 

since in this case the axinos inherit the thermally produced neutralino number density. 
The neutralino-to-axino decay offers a mechanism to shed large factors of relic density. For 
a case where m^o ~ 100 GeV and n^oh^ ~ 10 (as can occur in the mSUGRA model at 
large tuq values) an axino mass of less than 1 GeV reduces the DM abundance to below 
WMAP-measured levels. 

The lifetime for these decays has been calculated, and it is typically in the range of 
r(xi 07) ~ 0.01 — 1 sec[Q. The photon energy injection from Xi ~^ 07 decay into 
the cosmic soup occurs typically before BBN, thus avoiding the constraints that plague 
the case of a gravitino LSPll^. The axino DM arising from neutralino decay is generally 
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considered warm or even hot dark matter for cases with ~ 1 — 10 GeV[^]. Thus, in the 
mSUGRA scenario considered here, where ~ 1 — 10 GeV, we usually get warm axino 
DM from neutralino decay. 

2.3.3 Thermal production of axinos 

Even though axinos may not be in thermal equilibrium in the early universe, they can still 
be produced thermally via scattering and decay processes in the cosmic soup. The axino 
thermally produced (TP) relic abundance has been calculated in Ref. ||2^, and is given 



in Ref. |31] using hard thermal loop resummation as 



where gs is the strong coupling evaluated at Q = Tr and N is the model dependent color 
anomaly of the PQ symmetry, of order 1. For reference, we take gsiTR = 10^ GeV) = 0.932 
(as given by Isajet 2-loop gs evolution in mSUGRA), with gs at other values of Tr given 
by the 1-loop MSSM running value. The thermally produced axinos qualify as cold dark 
matter as long as ms ~ 0.1 MeV|29, 31|. 

In Fig. H, we plot the re-heat temperature needed to thermally produce various abun- 
dances of axinos versus the Peccei-Quinn scale fa/N. We plot values of Qj^h"^ = 0.001 
(solid), 0.01 (dashed) and 0.1 (dot-dashed), assuming values of nia = 10~^ (purple), 10~^ 
(green) and 1 GeV (maroon). We only plot solutions with Tji ~ 10^ GeV, since for lower 
values of Tr ~ 10^"'^ GeV, Eq. is expected to break down. We see from the curves that 
in order to achieve T/j values ~ 10^ GeV, we will need values of fa/N on the large side: 
10^^ — 10^^ GeV. We also see that the purple curves- with lowest values of ma 100 
keV- give the largest Tr values. Of course, from the preceeding discussion, large values 
of fa/N also give more axion dark matter, independent of any other parameters. Thus, 
to achieve high values of Tr, we will likely need to examine scenarios with mostly axion 
CDM, combined with smaller amounts of cold and warm axinos. 

3. Preferred mSUGRA parameters with mainly axion CDM 

In this section, we generate sparticle mass spectra using the Isasugra subprogram of the 



event generator Isajet[16|. Isasugra performs an iterative solution of the MSSM two-loop 
RGEs, and includes an RG-improved one-loop effective potential evaluation at an optimized 
scale, which accounts for leading two- loop effects [ 32 1. Complete one-loop mass corrections 



for all sparticles and Higgs boson masses are included pSj.^ 

Our first results are shown in Fig. where we examine the mSUGRA point with 
(mo, mi/2) ^0) tan /?, sgn{ij)) = (1000, 300, 0, 10, -|-1) (where all mass parameters are in GeV 
units). We also take mt = 172.6 GeV. For this point, the neutralino relic density computed 



^The case of < m^o was recently examined in Ref. [^2|. In their results, they always take ^ah^ ~ 0. 
We have checked using the Micromegas program|l^] (to calculate the stau relic density, which is not handled 
by IsaReD) that the value of Tr generated in the stau NLSP region is always less than the corresponding 
values generated in the neutralino NLSP regions for the cases considered in this section. 
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nL^h- = 0.001 (solid), 0.01 (dashed), 0.1 (dashed-dotted) 
nij (GeV) = 0.0001 (purple), 0.01 (green), 1 (maroon) 




f^/N (GeW) 

Figure 1: A plot of the expected re-heat temperature of the universe versus PQ breaking scale 
fa/N for nj^h'^ = 0.001 (solid), 0.01 (dashed) and 0.1 (dot-dashed), and with ma = 10"^ (purple), 
10~^ (green) and 1 GeV (maroon). (The solid purple and dot-dashed green, and also the solid green 
and dot-dashed maroon lines coincide.) 



by IsaReD|34] is ^y.oh = 8.9, so the point would be excluded under the assumption that 
thermal neutralinos make up the dark matter. In frame a)., we plot the values of Oa/i^, 
ri|'^/i^ and Q.^^^h'^ versus fa/N under the assumption that Tr = 10^ (dashes), lO'^ (solid) 
and 10^ GeV (dot-dashed). We assume the axion relic density is as given by the central 
value of Eq. |2]|. We require as well that the sum Q.ah? + ^l^h? + Q^^^h^ = 0.11, i.e. 
that the combination of three components of axion and axino DM saturate the WMAP 
central value. For each value of fa/N, the value of needed to saturate the measured 
DM abundance is calculated, and listed in frame b). in GeV units, along with in eV 
units. The axion abundance is of course independent of Tr. The abundance of Vl^^h? is 
fixed mainly by requiring the total relic abundance saturate the measured central value, 
and since Tr is fixed, this means we can compute the needed value of m^. For low values 
of fa/N ~ 10^^ GeV, the DM abundance is dominated by Q'F^h'^. (the curves for ^'^^h'^ 
for all three cases of overlap to within the line resolution). But comparing with frame 
b)., we see for almost all of this range, rria < 100 keV, meaning the bulk of axino DM 
is actually warm, in contradiction to what is needed to generate large scale structure in 
the universe. An exception occurs in the case of Tr = 10^ GeV (barely enough for non- 
thermal leptogenesis) , where rria moves to values higher than 100 keV. At the highest 
fa/N - 3 X 10^^ GeV, axion CDM dominates the relic abundance. In this case, rria must 
drop precipitously so that Qj^h"^, which wants to rise with increasing fa/N, instead sharply 
drops. The region with mainly axion CDM is robust in that it gives rise to a consistent 
cosmology for all choices of Tr: for this case, the axino mass can drop below 100 keV into 
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Figure 2: Axion and TP and NTP axino contributions to dark matter density for Tr = 10^ GeV, 
10^ GeV and 10^ GeV versus PQ breaking scale fa/N. 



the warm DM region, since now axinos will only be a small component of the DM density. 



In Fig. ^, we plot again the same quantities as in Fig. ^, but this time we keep ma 
fixed to a value of 100 keV (solid) and 1 MeV (dashed), and we allow Tr to vary in order 
to maintain the WMAP measured abundance of CDM. Frame a), shows the relic density 
of all three components of axion/axino dark matter, while frame h). shows the value of 
Tr needed for each value of fa/^ ■ We see that for low values of fa/N, the value of Tr is 
weh below the 10*^ GeV regime, and in fact doesn't even exceed 10^ GeV for the case of 
rria = I MeV. In the case of rria = 100 keV, Tr exceeds 10^ GeV for fa/N ~ 10^^ GeV, 
and approaches a maximum for the case of mainly axion dark matter. 

Next, we explore the mSUGRA vs. mx/2 plane for the presence of solutions with 
Tr ^ 10^ GeV so they yield consistent baryogenesis mechanisms. In our first try, we set 
fa/N = 1.2 X 10^^ GeV so that the measured dark matter density is saturated by cold axions 
(we assume the factor of three downward fluctuation in l^a^'^i which allows for an increased 
value of fa/N). We will assume equal portions of TP and NTP axinos, which saturate the 
1-a error bars on the WMAP measured QcDuh^ value: ^^^h^ = Vt^'^^h'^ = 0.003. We 
also adopt mSUGRA parameters Aq = 0, tan /? = 10 and ^ > 0. Using these values, we 
calculate the sparticle mass spectrum, Q-oh? and m-o at each point in mSUGRA space. We 
then determine the necessary value of rria (from fi^o /i^), and then calculate the required 
value of Tr (from ^tF^h^). We plot in Fig. ^ the color-coded regions of Tr values, along 
with contours of logj^g ^R- The lower right red region is excluded due to lack of appropriate 
EWSB, while the left-side red region yields a stau NLSP. The gray region is excluded by 
LEP2 limits on the chargino mass. 
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Figure 3: Axion and TP and NTP axino contributions to dark matter density for ma = 100 keV 
versus PQ breaking scale fa/N . 



From Fig. |^ we see that the usual regions preferred for neutrahno cold dark matter 
actually give the lowest values of Tr: we find Tr < 10^ GeV in the stau co-annihilation 
region and in the hyperbolic branch/focus point (HB/FP) region. These values may even 
be too small to sustain electroweak baryogenesis. The central blue regions of the plot 
accommodate the largest values of Tr ~ 10^ GeV. While yielding a much higher Tr value 
than the stau and HB/FP regions, even these regions do not yield a high enough Tr value 
to sustain non-thermal leptogenesis. 

Fig. ^ shows the corresponding contours of rria- The values range from rria ~ 100 (600) 
GeV in the stau (HB/FP) regions to values of rria < 0.05 GeV in the regions of high Tr. 

To push the value of Tr higher, we would need to diminish even further the value 
of rua (as suggested by Fig. §6).) which also diminishes the amount of NTP axino dark 
matter. In Fig. |6|, we again plot the mo vs. mij2 plane for the same parameters as in Fig. 
^ so that we saturate the CDM abundance with axions. However, in this case we adopt 
Vll^h'^ = 0.006 (saturating the WMAP QcDAih^ error bar), and take Q^^^h^ = 6 x 10"^. 
The much smaller value of Q.^'^^h? (than that used in Fig. Q) means the value of will 
be much smaller all over the plane than in the Fig. ^ case. To balance the lower value of rria 
in the fixed value of Q,'^^h'^, a much higher value of Tr will be needed. We now see that the 
contours of Tr plotted in Fig. |^ move well into the 10^ GeV regime: enough to sustain the 
non-thermal leptogenesis mechanism. In fact, the preferred regions of high Tr are precisely 
those regions of mSUGRA parameter space that are most disfavored by neutralino CDM! In 
this case, the stau and HB/FP regions, preferred in mSUGRA, can only sustain Tr values 
in the 10^ GeV range. The regions with mo ~ 1 TeV and mi/2 ~ 200 — 400 GeV allow for 
Tr values well in excess of 10^ GeV. This region of parameter space is usually neglected in 
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(GeV) 

Figure 4: Contours of constant Tji in the mo vs. vfiiji plane for Aq ~ 0, tan/3 — 10 and /i > 0. 
We assume Vljr? = 0.11, and i^J^h^ = ^^'^^ = 0.003. The lower right red region is excluded due 
to lack of appropriate EWSB, while the left-side red region yields a stau NLSP. The gray region is 
excluded by LEP2 limits on the chargino mass. 



simulation studies for the LHC, since it severely disagrees with the conjecture of thermally 
produced neutralino CDM. For this reason, we list a benchmark point A in Table |l| with 
niQ = 1500 GeV, mi/2 = 200 GeV, = 0, tan/? = 10 and fi > 0. In this region, squarks 
and sleptons have mass in the TeV range, while nig ~ 500 GeV. LHC collider events should 
thus be dominated by gluino pair production, followed by gluino three-body decays into 
ggXi and qq'xf final states. The xf and X2 ^i^^ decay into ff'xi and ffxi respectively, 
where / denotes any SM fermion states whose decay modes are kinematically allowed. In 
particular, the decays X2 ~^ Xi^~^^~ and X2 ~^ Xif^~^IJ'~ each occur at ~ 3%, since the 
decay is dominated by Z* exchange. 

Fig. ^ shows contours of rria for the case as shown in Fig. |6|. Here, we see much 
smaller rria values below 100 keV are generated. These low values of nia would yield warm 
thermally produced axinos. However, since the bulk of DM is constituted by axions, the 
temperature of the small fraction of axinos is not relevant. 

In Fig. ^, we again plot contours of constant in the mo vs. mi/2 plane for ^0 = 
and n > 0, but this time for tan/3 = 30. The larger value of tan/3 leads to larger 
values of b and r Yukawa couplings, and a lower value of myi|35]. This in turn leads to 
larger rates for neutralino annihilation via s-channel A* exchange diagrams, and somewhat 
lower neutralino relic density values. We again asume = 0.11, ^l^h"^ = 0.006 and 

rj^^^ = 6 X 10^^ and calculate the requisite value of Tr. From the figure, it is seen that 
the stau and HB/FP regions lead to lower values of Tr ~ 10^ — 10^ GeV, while regions 
with niQ ~ 800 — 2000 GeV can still lead to Tr in excess of 10^ GeV, enough to sustain 
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Figure 5: Contours of constant in the toq vs. mi/2 plane for Aq = 0, tan/5 = 10 and /i > 0. 
We assume Oa^i^ = 0.11, and nl^K^ = f7f = 0.003. 




(GeV) 



Figure 6: Contours of constant Tr in the mo fs. mi/2 plane for Aq — 0, tan/? = 10 and /x > 0. 
We assume nah"^ = 0.11, and f^P/i^ = q.OOG and = 6 x 10"^. 

non-thermal leptogenesis. 

The corresponding contours of for the tan /3 = 30 case are shown in Fig. ^. They 
range from above 1 GeV in the stau co-annihilation and HB/FP region, to below 100 keV 
in the regions of Tr > 10^ GeV. We list in Table | an mSUGRA benchmark point B with 



- 11 - 



1500 



> 

dJ LOOO 



500 



0!0003^-f ^ 




m-^ (QeV) 

1>0.0013 



0.001 1 



1000 2000 3000 4000 5000 6000 



0.0009 



0.0007 



0.0005 



0.0003 



0.0001 



CGeV) 



Figure 7: Contours of constant ma in the mo vs. mi/2 plane for Aq = 0, tan/? = 10 and fi > 0. 
We assume Qah^ = 0.11, and flT^h^ = 0.006 and fl^^^ = 6 x 10"^. 
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Figure 8: Contours of constant Tn in the mo vs. mi/2 plane for Aq — 0, tan/? ~ 30 and /i > 0. 
We assume flah^ = 0.11, and flTPh^ ^ 0.006 and Q^^^ = 6 x IQ-*^. 



tan P = 30 and mainly axion CDM. In this case, as in the case of benchmark point A, 
gluino pair production occurs at a large rate. However, in this case, X2 ^ Xi^ at ~ 100%, 
so LHC events will be rich in multi-jet plus Z plus EJp^^^ signatures! 36 1. 

It should now be apparent that by conjecturing a large value of the PQ scale fa/N 
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Figure 9: Contours of constant in the uiq vs. mi/2 plane for Aq 
We assume Qah^ = 0.11, and flTPh^ = 0.006 and f7f = 6 x 10"^. 



0, tan/3 = 30 and fj, > 0. 



such that axion CDM saturates the rehc density, then by taking decreasingly low values 
of axino mass, large values of Tr may be generated. To illustrate this graphically, we plot 
in Fig. IC the value of Tr required versus rria for the mSUGRA point mo = 1000 GeV, 
300 GeV, ^0 = 0, tan/3 = 10 and /i > 0. We assume Qah'^ = 0.11, and that 
0.006. For nid = 0.08 GeV, the axino portion of the relic density is 



1^1/2 = 



NTP 



comprised almost entirely of non-thermally produced axino DM from neutralino decay. As 
we decrease from this value, the portion of NTP axino DM decreases, and since the 
sum is constant, the TP portion increases. Since the value of Q.'F^h'^ is proportional to rria, 
a large increase in Tr is needed to keep pace. We see in the extreme limit, values of Tr as 
high as lO^'^ GeV can be generated (putting us in conflict with the gravitino problem) for 
values of rria as low as 10~^ GeV (far below the usual theory expectations for ma[29|). 



4. Summary and conclusions 

In this paper, we have examined the consequences for the mSUGRA model if dark matter 
is composed of an axion/axino admixture, rather than neutralinos. We have considered 
this scenario along with cosmological consequences of the gravitino problem (which re- 
stricts 771-3/2 > 5 TeV and re-heat temperatures Tr ~ 10^ GeV) and leptogenesis. While 
thermal leptogenesis requires Tr ^ 10^ GeV (in conflict with the gravitino problem), non- 
thermal leptogenesis- wherein heavy right-hand neutrino states are produced additionally 
via inflaton decay- can allow for successful baryogenesis with Tr ~ 10^ GeV. (In addition, 
Affleck-Dine leptogenesis may occur at these values of Tr, although it may also occur at 
even lower Tr values.) 
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parameter 


Pt. A 


Pt. B 


mo 


1500 


1000 


mi/2 


200 


300 


Ao 








tan P 


10 


30 


fJ' 


304.5 


368.4 


m~g 


OOo.Z 


1 1 0.0 




1541.1 


1178.4 


Til 7 

"Hi 


912.3 


774.2 




1 ORA ^ 
iZD4.0 


QRK 7 




ibUU.D 


iUUo.o 




1 /I Q 7 
145. ( 


007 Q 




148.0 


227.0 


Al 


80.0 


122.4 




1510.6 


912.2 




112.4 


112.8 


BF{b 57) 


9.2 


6.5 


3.1 X 10-"^ 


2.5 X 10- 



Aafy^^ 1.5 X 10-1° 9.7 x 10"^° 

Table 1: Masses in GeV units and parameters for two mSUGRA model benchmark points with 
mainly axion CDM: l^^/i^ = 0.11. We take fa/N = 5 x 10" GeV, along with = 0.006 and 
^NTP ^ g X iQ-6 j^^gQ ^ ;^72.6 GeV. 



We explored mSUGRA parameter space for regions of high Tr with three components 
of dark matter: dominant axion CDM, along with small portions of thermally and non- 
thermally produced axino DM (which may be either warm or cold). We find the highest 
values of Tr occur in the regions of mSUGRA space which are typically most disfavored 
by neutralino CDM. Likewise, the regions of mSUGRA parameter space most favored by 
neutralino CDM are actually most disfavored by mixed axion/axino DM. By combining 
high Tr values with fine-tuning considerations (which prefer lower values of rriQ and espe- 
cially ?Tii/2)j we find mSUGRA with mainly axion CDM prefers uiq ~ 800 — 2000 GeV, 
with mi/2 ^ 150 — 400 GeV. There are several consequences of this scenario: 

• LHC SUSY events will be dominated by gluino pair production, followed by gluino 
three body decays to charginos and neutralinos. 

• Current and future WIMP direct and indirect dark matter detection experiments will 
likely find null results. 



The ADMX|37], or other direct axion detection experiments, stand a good chance 
of finding an axion signal. The ultimate axion rate predictions are of course model 
dependent. 
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Figure 10: Plot of the value of needed versus mg for the mSUGRA point with mg = 1000 
GeV, mi/2 = 300 GeV, Aq = 0, tan/? 10 and ^ > 0. We assume ^lah? = 0.11, and that 
+ n^'^P = 0.006. 



A related scenario for sparticle spectra with mainly axion CDM has been put forward 
in Refs. |^8| in the context of t — 6 — r Yukawa-unified SUSY models. These models- 
expected from simple 50 (10) SUSY GUT theories- predict a spectrum of scalars in the 
range of 5 — 15 TeV: much higher than allowed in mSUGRA. Additionally, Yukawa unified 
models require tan /? ~ 50. Thus, the mSUGRA model- in contrast to Yukawa-unified 



models|39]- allows for the possibility of first and second generation squark masses which 
are accessible to LHC, and which would augment the gluino production rate, since q ^ qg 
decay is expected. In addition, mSUGRA models with mainly axion CDM allow for lower 
values of tan /? and hence smaller values of b and r Yukawa couplings. In this case, a lower 
multiplicity of 6-quark jets is expected in LHC SUSY events. 
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